Equation of State of Oscillating Brans-Dicke Scalar and Extra Dimensions 
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We consider a Brans-Dicke scalar field stabilized by a general power law potential with power 
index n at a finite equilibrium value. Redshifting matter induces oscillations of the scalar field 
around its equilibrium due to the scalar field coupling to the trace of the energy momentum tensor. 
If the stabilizing potential is sufficiently steep these high frequency oscillations are consistent with 
observational and experimental constraints for arbitrary value of the Brans-Dicke parameter u). We 
study analytically and numerically the equation of state of these high frequency oscillations in terms 
of the parameters uj and n and find the corresponding evolution of the universe scale factor. We 
find that the equation of state parameter can be negative and less than -1 but it is not related to 
the evolution of the scale factor in the usual way. Nevertheless, accelerating expansion is found for 
a certain parameter range. Our analysis applies also to oscillations of the size of extra dimensions 
(the radion field) around an equilibrium value. This duality between self-coupled Brans-Dicke and 
radion dynamics is applicable for w = — 1 + 1/D where D is the number of extra dimensions. 
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I. INTRODUCTION 

Observations of the magnitude-redshift relation of dis- 
tant Type la supernovaejj and CMB measurements 
have indicated that the universe is undergoing a period 
of accelerated expansion. Attributing this acceleration 
to a cosmological constant is not completely satisfac- 
tory (see e.g. Q) because fine tuned initial conditions 
are required in order to solve the coincidence problem 
(why the vacuum energy is dominating the energy den- 
sity right now). Thus, acceleration is attributed to the 
influence of a redshifting non-luminus form of energy 
with negative pressure called by many authors 'dark en- 
ergy' or 'quintessence' The origin of this form of en- 
ergy remains unknown. Several types of scalar fields |4| 
have been proposed as sources of dark energy includ- 
ing the dilaton, the inflaton, supersymmetric partners 
of fermions, Brans-Dicke (hereafter BD) scalarsQ, etc. 

The dynamical evolution of these scalars is determined 
by specially designed potentials so that the scalar field 
energy density tracks the radiation-matter component at 
early times while at late times it dominates and leads to 
acceleration. The main drawback of these models is that 
there is usually no physical motivation for the proposed 
scalar fields and potentials. 

Cosmological theories with extra dimensions generi- 
cally contain a scalar field of geometrical origin which 
describes the size of extra dimensions. This is known 
as the radion field. These theories have the potential of 
providing a physically motivated solution to the hierar- 
chy problem by postulating that the fundamental Planck 
mass M* is close to the TeV scale 0, 0,0 . This is possi- 
ble in theories with extra dimensions Q because Gauss's 
law relates the Planck scales of the 4 + D-dimcnsional 
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where 60 is the present stabilized size of the extra dimen- 
sions. 

According to observations, the internal space should be 
static or nearly static at least from the time of primor- 
dial nucleosynthesis 0- This means that at the present 
evolutionary stage of the Universe only small fluctua- 
tions over stabilized or slowly varying compactification 
scales (conformal scales/geometrical moduli) are possi- 
ble. Thus the size of the extra dimensions b(t) (the radion 
field) is usually assumed to be stabilized to its present 
value bo by a 'radion stabilizing potential' V(b). Obser- 
vational constrains require this stabilization and theoret- 
ical mo dels |1 Of have been proposed to justify it. Within 
the framework of multidimensional cosmological models 
the role of the size of extra dimensions as a stabilized 
scalar field with possible excited states was investigated 
in [Til ll^ . Il3| where they were called gravitational exci- 
tons. 

In a recent paper we demonstrated that a generic 
feature of these theories is the existence of radion oscilla- 
tions (induced by redshifting matter density) around the 
equilibrium value feo at the minimum of the stabilizing 
potential. It was also shown that these oscillations can 
lead to cosmological periods of accelerating expansion of 
the universe. However the radion mass required for such 
extremely low frequency oscillations was too small (at 
the classical level) to be consistent with fifth force exper- 
imental constraints [TEl In the same paper^ij it was 
also demonstrated that the dynamics of the radion field 
are formally equivalent with the dynamics of a massive 
BD scalar with u> = — 1 + 1 AD where D is the number of 
extra dimensions (see also [rH ITc| ). 

Quintessence models based on BD scalars (extended 
quintessence) have the theoretically appealing property 
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that the same scalar field that participates in the gravity 
sector is enhanced by a potential 19] to play the role of 
quintessence. Even though no oscillating BD scalars have 
been studied so far, there is an extensive literature on 
extended quintessence |20j, |2jJ, |22( . 

Here we focus on the cosmological evolution of high 
frequency oscillations ^| 23] of a BD scalar^ (or equiva- 
lently a radion field) stabilized by an arbitrary power law 

potential. For large enough frequencies v ~ 10 12 Hz these 
oscillations are consistent with fifth force constraints fl6l 
l23l ] and they can play the role of dark matter and/or dark 
energy for a wide range of parameters. 

The structure of this paper is the following: In sec- 
tion II we review the cosmological evolution of an os- 
cillating minimally coupled scalar field. In section III 
we demonstrate the equivalence between radion dynam- 
ics in D extra dimensions and massive BD scalar <fr with 
uj = — 1 + l/D, 4> — b 2D and focus on the dynamics of 
an oscillating BD scalar. We derive a virial theorem that 
connects the kinetic and potential energies of the oscil- 
lating field in terms of the parameters u> and n. This 
theorem is verified by numerically solving the dynamical 
cosmological field equations for various parameter values. 
We then use this virial theorem to find the equation of 
state that connects the effective pressure with the effec- 
tive energy density of the oscillating field. In section IV 
we use the equation of state to find the redshift rate of 
the effective energy density and the evolution of the cos- 
mological scale factor. We find that the energy redshift 
and the scale factor evolution are independent of lj and 
depend only on the stabilizing potential power index n. 
The validity of these results is demonstrated numerically. 
Finally in section V we conclude and discuss possible ex- 
tensions and open issues related to this work. In what 
follows we set 8nG = 1. 



II. MINIMALLY COUPLED OSCILLATING 
SCALAR 

As a warm up exercise we evaluate the equation of state 
for a minimally coupled coherently oscillating scalar in an 
expanding background. This analysis was originally done 
in Ref. [24J but we briefly review it here for completeness 
and for comparison of the results with the BD scalar 
oscillations analysis performed in the following sections. 

Consider a minimally coupled scalar field whose dy- 
namics is determined by the Lagrangian density 



c = -d^d^-v^) 



(2.1) 



Assuming a Robertson- Walker metric we obtain the field 
equation 



3^ = -^, 



(2.2) 



obtain the energy density p and the pressure p as 

p = \i> 2 + v{4>) 
p = \fi - V{4>) 

Using equations l|2.2[) and l|2.3[) we obtain 



dt 



p = -3Hc/) 2 



We now consider a potential of the form 
V(<£) = A|C 



(2.3) 
(2.4) 

(2.5) 
(2.6) 



where n > and focus on oscillations of the scalar field 
around its minimum. We also define the mean value of a 
time dependent quantity q(t) as 



(2.7) 



where T is the period of one oscillation. The mean equa- 
tion of state of the oscillating field is 



p = wp 



(2.8) 



We now evaluate the equation of state parameter w in 
terms of n. 

Using equations (|2.3() , (|2.4H , 12.5(1 and l|2.8|l we obtain 

^ = -3H4> 2 = -3# 7 p = -3-7p (2.9) 
dt a 

where 7 = w + 1 and we have assumed T <C H~ l . Thus, 
on timescales At <C H^ 1 we have ^ -C 1 and we can 
set 

P=\<t> 2 + V{4>) ^ V{(p m a X ) = Vmax = P (2.10) 

where <fimax is the amplitude of the 4> oscillations and 
V max is the corresponding maximum potential energy. 
We now have 



p + p 



V2 



p P T„ 



o V, 
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We also use equation l|2.10|l to find the time T„ 
it takes the field to go from to 4> m ax 



(2.11) 



T 



d<p{l-V/V max )- 1 / 2 (2.12) 



Using equations l|2.11|l and (|2.12|) we find 

J*-* #(1 - U/Vw) 1 ^ 2r 



7 = 2- 



where H 



From the energy momentum tensor we 



#(1 - V/V max )-i/ 2 ™ + 2 
where we made use of (|2.6|) . 



(2.13) 
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The energy density redshift is now obtained from equa- 
tion 12. 9f) as 



P : 



a 

Pa \ — 

ao 
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Using this equation and the Friedman equation 
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we obtain 



t 



2/37 = ^ 



(2.14) 

(2.15) 
(2.16) 



and accelerated expansion is obtained in the range < 
n < 1. The cosmological effects of this accelerated ex- 
pansion for an oscillating minimally coupled scalar have 
been discussed in Ref. 25]. The scalar field evolution 
for power index n in the above range is well defined even 
though for n < 1, V'(4>) is weakly singular at <fi = 0. V(4>) 
can be made mathematically more appealing by the field 
redefinition 4> — > ((f) 2 + cb 2 )^ 2 , </> c — > which does not 
change the above results 25] . In the next section we will 
show that even though the equation of state is different 
for an oscillating BD scalar, accelerated expansion is ob- 
tained for the same range of the power index n. 



III. VIRIAL THEOREM AND EQUATION OF 
STATE FOR OSCILLATING RADION AND 
BRANS-DICKE SCALAR 

In this section we demonstrate the equivalence between 
BD and Radion dynamics (for specific values of the BD 
parameter us) and then focus on the BD theory to derive 
a virial theorem and the equation of state of the scalar 
field oscillations. 

Consider the BD action 



dW^S <j>R--d^d> 1 <t>-2V{4>) +L flaid 

(3.1) 

where Lfi U id describes matter and/or radiation. We as- 
sume a flat (k = 0) Friedmann-Robertson- Walker model 
whose metric reads 



ds 2 = -dt 2 +a 2 (t)[dr 



l dd 2 



r sm 



Bd<\> 2 \ (3.2) 



The combination of equations i|3.1|) and (|3.2() leads to the 
dynamical equations for and a 



a 4> 
a 4> 
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(3.3) 



(3.4) 



(3.5) 



From the field equations we can read the effective energy 
and pressure (see e.g. |20ll26j ]') for the field which end up 
being, 



U 2~6 +y -^ 
Am a 



and 



V + (j) + 2-4> 
a 



(3.6) 



(3.7) 



The BD theory is interesting in its own right as a gener- 
alization of Einstein gravity. In addition however many 
physically motivated theories reduce to BD gravity for 
specific values of ui. For example the superstring dilaton 
is akin to a BD theory with u> = — 1. Also Einstein grav- 
ity generalized to 4+D dimensions reduces to an effective 
4-dimcnsional BD theory with u> = — 1 + 1/D. This can 
be demonstrated as follows: 

We consider a 3-brane embedded in a toroidally com- 
pact (4 + D) dimensional space-time R 1 x R 3 x T D with 
the D extra spatial dimensions (the bulk) stabilized at a 
size bo. The total action may be written as 



Stot — S\y U lk + Sbrane 



(3.8) 



where 

Sbulk 

and 



d i+D xV-g(*+n) {-±=n[g\ + C hulk ) (3.9) 
\lD7rG / 



^brane 



(3.10) 



In 13.9fl £buik is the Lagrangian of the bulk fields apart 
from the graviton. These fields give rise to the stabilizing 
potential V discussed below. The background metric for 
the a toroidally compact (4 + D) dimensional space-time 
which is consistent with the symmetries of the brane-bulk 
system may be written as 



9mn = diag[l,-a 2 (t)g i:j ,-b 2 (t)g mn \ 



(3.11) 



where M, N run from to D + 3; i,j run from 1 to 3 and 
m, n run from 4 to D + 3. Also a(t) is the scale of the 
non-compact 3-dimensional flat space (the scale factor of 
the universe) and b(t) is the radius of the compactified 
toroidal space (the radion field). 

We assume that matter on the brane (Sbrane) is repre- 
sented by a classically conserved perfect fluid energy mo- 
\ v = (p+P^^Uu+pg^u (m> f = 0, ...,3) 



mentum tensor T, 
with V^T^ 



where it M is a future-oriented time-like 
vector — (l,03 + £>) in the basis H3.11JI and p is the 
energy density of the brane matter. Also p — wp is the 
corresponding pressure. We also include a radion stabi- 
lizing potential 



V(b) = b(t) D / d 4+D x C 



-bulk 



(3.12) 
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induced by some non-specified bulk dynamics £t>uikE3- 
For example a uniform bulk cosmological constant is rep- 
resented by V = b D A while a Casimir type potential 
would be V(b) = kb' 4 H3- We assume that bulk dy- 
namics lead to a potential V(b) that stabilizes b at 

/ M 2 \ 1/D 



with a vanishing cosmological constant. The equations of 
motion for the coupled (a(i), b(t)) system can be written 
as 28] 



b 

b ' 

a 
a 



2 S 



^ h 1 

'ab~W 



2V 



D 



a b 



-2 
dV 



D(D 
D 



2{D 



2) db 
I — 



2(D 



dV 



2) db 
—rrV H 



2(D + 2) 

P+(£>-l)p 
2(£» + 2) 



(3.14) 



In the picture of Refs. (see also |2g) (ADD) and 
(RS), the matter-radiation energy density is assumed to 
be localized on the brane corresponding to the a(t) scale 
factor. This localized energy density in general distorts 
the geometry of the compactified /^-dimensional space 
(the bulk), but as far as the overall properties and the 
evolution of the radion are concerned, it is correct to 
treat the energy density on the wall as just being aver- 
aged over the whole space as done on the RHS of these 
equations. This assumption is consistent with the results 
of references [2|| [3(J where it was shown that the cou- 
pling of the radion field to the energy momentum tensor 
is given generically through the trace (p — 3p) plus terms 
involving the stabilizing potential V(b). 

Equations (|3.14|) have been analyzed in the context of 
inflation (away from the stabilization point bo) in j2^ 
(ADD model). 

Similar equations j3Cllarise in the context of Randall- 
Sundrum (RS) models jg, where the hierarchy problem 
is solved using extra dimensions of smaller sizes at the ex- 
pense of introducing a non-flat background metric along 
the extra coordinates and a pair of branes whose distance 
b(t) is stabilized at bo by the potential V(b). 

From equations l|3.14|l it is clear that radiation (p ra d = 
\prad) has no effect on the dynamics of the radion. This 
is not the case however for matter. Redshifting mat- 
ter plays the role of a driving force and can induce ra- 
dion oscillations during both the matter and radiation 
eras. These oscillations backreact on the scale factor a(t) 
through the effective Friedman equations 1(3. 14|) and can 
affect the expansion rate. 

It is interesting to compare the form of the cosmo- 
logical equations (|3.14|) with the non-conventional cos- 
mology on the brane of Ref. [3l| which has shown that 
terms of 0(p 2 ) appear on the Friedman equation due to 
the extra dimensions. These terms can also be obtained 
from the system (I3.14f) in the limit (|) 2 <C p. In that 



limit, it is a good approximation to set b = in the 
second of equations 1(3.14(1 . For a stabilizing potential 
of the form V ~ (b — bo) n we thus obtain V ~ p™- 1 
and for a standard parabolic potential (n = 2) we ob- 
tain an extra p 2 contribution due to the potential in the 
first equation of the system (|3.14|) similar to the p 2 terms 
of Ref. |3jJ. For a generalized stabilizing potential the 
extra term in the Friedman equation is of the form p a 
(a = ^zj) and can provide a physical motivation for 
Cardassian expansion |32j. In what follows however we 
will not use the above approximation which ignores the 
effects of radion oscillations and we will focus on the cos- 
mological expansion induced by these oscillations for a 
generalized stabilizing potential. 

It is straightforward to show that we can obtain the set 
of dynamical equations l|3.14|l for the radion by setting to 
the corresponding BD set of equations (|3.3|l . H3.4(l . 1)3.5(1 

J. - l-D 



UJ 



■'4 



Thus the radion dynamics is equivalent to BD dynam- 
ics for particular values of u. Even though our study 
will be based on the general BD field equations l(3.3|l . 
(|3.4|) and (|3.5|l the main physical application of our re- 
sults will be A + D dimensional Einstein gravity for which 
uj = —1 + 1/D. Even though these values of u> are of 
O(l) they are consistent with experimental and observa- 
tional constraints because of the presence of the poten- 
tial V. This stabilizing potential can confine the effects 
of any modifications of Newton's law to scales smaller 
than 1mm 0, 0, |2^| . On these scales there are no firm 
experimental constraints on the form of the gravitational 
potential |l5j. 

We now consider the BD field equations (|3.3|) . (|3.4|l . 
(|3.5|) with a potential V((f>) having a minimum at <po- It 
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is clear from equation l|3.5|l that even if we consider an 
initially static field <f> = (f>o the redshifting matter will 
induce small oscillations around the equilibrium position 
(approximately <f>o) by slowly shifting the equilibrium po- 
sition of the field (f>. We will thus study the equation of 
state of the energy density of these small oscillations as- 
suming a potential of the form 



v^) = Ai0-0 O r = Ai<50r 



(3.15) 



Since not all three equations of the system l|3.3[) . I|3.4|) and 
(|3.5|) are independent (because of the Bianchi identities) 
we focus on (I3.3|) and l|3.5|) and write them in dimension- 
less form as 



1 , 1 



a? a 4> uj (j) 2 .. ... 
a z a (f> 6 <p z <p a A 



3 i 1 

a (2w + 3) [ a 3 ~ 



(3.16) 



(3-17) 

where we have assumed the presence of matter redshifting 
like 



Pm(t) = PQ, 



(—) 3 

a{t) 



(3.18) 



(radiation has no effect on the evolution of since p. 
3p r — 0). We have also set a = a(t p ) 
the present time) and 



1 where t p is 



POn 



■t 



POm 



t 



Pv 



(3.19) 
(3.20) 
(3.21) 



The — sign (+ sign) in equation (|3.17|l is valid when <fi > 1 
(4> < 1). As in the case of the minimally coupled scalar of 
section II (f> is periodic on timescales At << H~ l , With 
no loss of generality we may therefore focus on the half 



period 4r for which 



1 > 0. We also assume 



8(f) <C 1. Thus equation l|3.17[) is written as 



^ + 2^T3 ( 
This implies that 

d ,8(j) 2 6p~ v 
dl^~2~ + 2lu + 3 



1 



2w + 3 a 3 



3^ 
a 



--gSift) = -3-1 

2w + 3a 3 ; a 



(3.22) 



(3.23) 



Therefore the effective potential that determines the dy- 
namics of the oscillating BD scalar is 

6 Pv x±n 3 1 xj. 



v eff {4>) = 



2uj 



2uj + 3 a 3 



(3.24) 



For n > 1 the effective potential is not minimized at 
6cj) = but at the time-dependent value 



Oeql 



1 



6np v a" 



(3.25) 



If the amplitude of oscillations S<p max is much larger than 
the shift of equilibrium point due to matter redshift then 
we may ignore the shift and assume that the oscillations 
take place around S<f> = 0. In the opposite case when 
S(f>eqi 3> S(f>max we may ignore the oscillations and assume 
that S(f> ~ 5<p eq i at all times. In the later case the energy 
density of the BD scalar is dominated by potential energy 
V which redshifts like 



V&) 



-3n/(n-l) 



(3.26) 



and is therefore subdominant compared to matter at late 
times. Thus in what follows we focus on the former case 
assuming 



> 



(3.27) 



It will be shown that for n < 4 this condition is an at- 
tractor in the sense that even if it is not realized at early 
times it becomes true at late times i.e. 5(j) max redshifts 
slower than Stfi eq i. Thus since the oscillations of 8<fi take 
place symmetrically around 8<f> = and since the period 
and amplitude are approximately constant on timescales 
much less than H^ 1 the mean value of 8(f> (and of 8(f)) 
vanishes and equation (|3.23|) implies 



8(f> 2 6p~v . 
2 2uj + 3 ( 



6p v 
2w + 3 



Proceeding in similar way as in the previous section 
(equations (|2.9|) - (|2.13|l ) we have 

p5<j, max JgJ ) ( 1 _ yeff /yeff} 1/2 9 
XJ.2 = oyeff Jo V / V rnax) = ye ff 

9 max jj*- dS(P(i - veff/v„tt)-^ max 

(3.29) 

Using now this result and equations (|3.24|) , (|3.28|l we find 



n + 2 



and 



yeff = y<^ff (i 

v * max V 



80 2 



. ' til ft -Y' 



(in 



yeff 



Pv8<ff' 



2uj + 3 



(3.30) 



(3.31) 



This result is a virial theorem for BD scalar oscillations 
connecting the mean kinetic and potential energies. We 
have numerically confirmed the validity of this theorem 
by solving the full system l|3.16|l . I|3.17|l for p~ v ^S> 1 (high 
frequency oscillations) from = -^i== to to = 1< We used 
initial conditions 



5(p(U) = 

8(p(ti) = S(f> m ax < 1 

a{U) - (2/3)- 2 / 3 t- 2/3 



(3.32) 
(3.33) 
(3.34) 



In Figure 1 we show the virial ratio _ f ,„ as a func- 
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FIG. 1: The virial ratio _ %,„ as a function of n for three 

Pv<>4> 

values of u) (u) — lower line, u) — —0.5, and u> = — 1 upper 
line). Notice that the virial theorem (straight lines) is violated 
for n ~ 4 




0.2 0.4 0.6 OS I 



FIG. 2: The oscillations of 8(j> do not even overlap with 5(j> = 
for n — 6. The equilibrium value 5(j> min is larger than 
the oscillation amplitude S<f> max and redshifts like a -3 ^ 5 ~ 
t~ 2 ^ 5 (matter energy dominated in the evolution). Thus the 
assumptions of the virial theorem are violated. 



tion of n for three values of uj. The points correspond 
to numerically obtained values while the lines are the 
corresponding analytical results (|3.31(l . For n < 4 the 
analytically obtained results are well verified by the nu- 
merics. For n > 4 however the agreement is not good 
and the virial ratio rapidly drops instead of increasing. 
This domination of the potential energy is due to the vi- 
olation of the condition i|3.27[l . The evolution of <fi for 
n = 6 shown in Figure 2 indicates that the oscillations 
of 5(j) do not even overlap with S(f> ~ and therefore the 
assumption 5<j) — used in deriving the virial theorem 
is not realized for n > 4. This may be understood by 
comparing the rcdshift rate of 

6<p min ~ a- 3 ^- 1 ) (3.35) 

(cf equation H3.25JI ) with that of 

Wma X - a- 6/in+2) (3-36) 

which will be derived in the next section (see equation 
(I4.3|l 1. It is easy to see that for n < 4, Scj) max eventually 
dominates while for n > 4 it is 5<j) min that dominates 
leading to violation of the virial theorem and vanishing of 
the virial ratio since the potential energy dominates. The 
fact that the oscillations of <p(t) are effectively around 4> = 
1 for n < 4 is demonstrated in Figure 3 which shows the 
evolution of (f>(t) for n — 1.5. The oscillation amplitude 
&4> max redshifts like a -6 /(" +2 ) as predicted analytically 
(see equation (|4.3|l ). The period of the oscillations varies 
slowly with time and it can be calculated using equation 
(|3.28l) in a similar way as for a minimally coupled scalar 
(equation (|2.12|) ). The result is that the oscillation period 
T redshifts like 

3(n-2) 

T~a^+^ (3.37) 

i.e. it decreases for n < 2 while it increases for n > 2. 
This result is confirmed numerically in Figure 4 where 
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FIG. 3: For n = 1.5 the BD scalar oscillations are approx- 
imately symmetric around <j> ~ 1 (lu = was used in the 
evolution but the result was found to be independent of ui) 

we show ^( ^1) ) vs ln(a). The numerically obtained 
points are in good agreement with the analytically pre- 
dicted corresponding lines for n = 1.5 (decreasing T), 
n = 2 (constant T) and n — 3.0 (increasing T). 

We are now in position to calculate the equation of 
state that connects the mean effective pressure with the 
mean effective energy density of the BD oscillating field. 
Using the rescaling of equations l|3.19(l - l|3.21(l and equa- 
tions ipO)| . (j3~7|i we find that for small approximately 
symmetric oscillations 

p ~ ^S4> 2 + frSfr (3.38) 
6 

p ~ ^Stf - p v 6j> n (3.39) 
Thus keeping the notation of section II we have 

= P*+P4, = |* (340) 
P4> P4> 
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FIG. 4: The oscillation period increases (decreases) for n > 2 
(n < 2). 



Using now equations (|3.29() - l|3.31|l we find 

2cu n 



7 : 



um + 2o> + 3 



u> + 1 



(3.41) 



FIG. 5: The equation of state parameter w is less than — 1 
for negative u>. 



w — 



ujn — (2iv + 3) 
ojn + (2uj + 3) 



(3.42) 



For a minimally coupled scalar this result would imme- 
diately imply that the redshift of the energy density is 
of the form of equation lj2.14|l with 7 given by Ij3.41|l . 
If that were the case, superaccelerati on pol l2l) {H > 
which may be favored observationally 33J) would be eas- 
ily obtained for negative values of to which lead to w < — 1 
(see Figure 5). For a BD scalar however this is not the 
case. It is easily seen by using e.g. equation (|3.23|) that 
the continuity equation 



■3-000- 
a 



-P4>) 







(3.43) 



is not realized I2C 



redshift like 



-3 7 



for p0 and p^,. Therefore does not 
The redshift rate of the effective en- 



ergy density and the corresponding evolution of the scale 
factor a{t) are studied in the next section. 



IV. 



ENERGY REDSHIFT AND SCALE FACTOR 
EVOLUTION 



Even though the continuity equation (|3.43|) is not re- 
alized for a BD scalar, equation l|3.23|l can be used to 
derive the effective energy density redshift. For symmet- 
ric oscillations this equation implies 



6p v 



-5<t) r > 



dV 2 ' 2w + 3~ a 
Using now equations (|3.28|) and l|3.29|l we find 



dt max 



. 6 £ V e//_!L_ 

max m 1 o 

a n + 2 



(4.1) 



(4.2) 



which implies 



V. 



eff 



and leads to equation 
Using now equations 
find 



* ~ 5tf nax (4.3) 

used in the previous section. 
. (|33UI) and 



we 



V eff ^ a ~^ 

r rn.n.T ^ 



(4.4) 



which gives the redshift rate of the BD energy density. 
It is interesting to note that even though the equation 
of state parameter w differs from that of a minimally 
coupled scalar and depends on u, the redshift rate of 
the energy density is independent of u> and is similar to 
the corresponding result of a minimally coupled scalar. 
This result is numerically verified in Figure 6 where we 
show a plot of —ln(fl<f,(a) / p<f,(a = 1)) vs ln(a) obtained 
numerically for various values of the scale factor a and for 
three values ofn (n — 1.5, n — 2.0 and n — 3). According 
to equation 14.41) , the numerically obtained points should 
lie on straight lines with slope 3r where 



r = 



2n 



(4.5) 



This is indeed realized as shown in Figure 6 where we 
also show these straight lines. 

We now focus on the evolution of the scale factor in 
the presence of an oscillating BD scalar given the above 
derived energy density redshift. To isolate the effects 
of the scalar we ignore the presence of matter at this 
stage. Thus the rescaled Friedman equation l|3.16|) is 
approximated by 



PvC 



Ba 



-3T 



(4.6) 
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FIG. 6: The energy density of the oscillating BD scalar red- 
shifts according to the analytical predictions for three values 
of n. 
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FIG. 8: The numerically obtained evolution of the 
in good agreement with the analytical prediction < 
f^for n = 1.5. 



dependent redshift of the minimum S(j) eq i of the effective 
potential (|3~24"|) . 

The scale factor evolution derived above when matter 
is not present can be combined with equation H3.36(l to 
find the time evolution of 54> max . The predicted power 
law Sipmax ~ t 2 / n is verified numerically in Figure 8 for 
n = 1.5. 



RADION OSCILLATIONS - OUTLOOK 



FIG. 7: The numerically obtained evolution of the scale factor 
is in good agreement with the analytical prediction for n — 1.5 
when matter is not present. 



where B is a constant and T determines the redshift rate 
of and is given by equation l|4.5l) . This implies that 



(4.7) 



Thus even though 7 7^ V in the case of a BD scalar os- 
cillations, the induced scale factor evolution is identical 
to that of a minimally coupled scalar shown in equa- 
tion H2.16fl . This result is numerically verified in Figure 
7 where we show a plot of ln(^jp-) vs In(j-) (we chose 
to ~ 2/3 and clq = a(ta)). The dark continous line corre- 
sponds to the numerical scale factor evolution for n = 1.5 
while the light dashed line is the corresponding analyti- 
cal result of equation l|4.7|l . The agreement is quite good, 
confirming the analytical expectation. It is interesting to 
note that this result implies accelerated expansion for 
< n < 1 as in the case of a minimally coupled scalar. 
The oscillations in this case however can be generically 
induced not just by initial conditions but also by the time 



The above derived results can easily be applied to the 
case of radion dynamics. By setting u> = — 1 + 1/D in 
equation i|3.42f> we can obtain the equation of state pa- 
rameter for high frequency radion oscillations in terms of 
the number of extra dimensions D as 



n(D - 1) + (D + 2) 
n(D- 1)- (D + 2) 



(5.1) 



Clearly w is negative and in fact less than — 1 for a wide 
range of parameters n, D. The evolution of the scale fac- 
tor however is independent of D and is given b y eq uation 
(|4.7|) . This result is in agreement with Ref. ljj where 
the case n = 2 was studied. 

Even though radion oscillations can not lead to super- 
acceleration, they do lead to accelerating expansion for 
n < 1 as in the case of a minimally coupled scalar field. 
Thus, the mechanism of Ref. pfij of quintessence from 
an oscillating minimally coupled scalar with n < 1 can 
be extended to the case of an oscillating BD scalar and 
in particular the case of an oscillating radion. In fact an 
oscillating radion has two potential advantages over an 
oscillating minimally coupled scalar. First its origin is 
better defined and better motivated since it is a generic 
geometrical prediction of theories with extra dimensions. 
Second, as discussed below, tracking can in principle be 
achieved at early times by the coupling to the matter 
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energy momentum tensor which generically induces ra- 
dion oscillations due to the matter redshift (the effective 
potential is time dependent). Since this possibility does 
not exist for a minimally coupled scalar, an exponential 
potential has to be used 25] for large (j> so that the energy 
density in <p tracks the matter component at early times 
while at late times the field oscillates around its poten- 
tial minimum with n < 1 thus producing the observed 
accelerated evolution of the scale factor. The potential 
used in Ref. [2{| that combines an exponential behavior 
at large <f> with power law behavior with power index n 
around its minimum is 

V(<f>) = V (cosh \4> - l)"/ 2 (5.2) 

Our results on the scale factor evolution indicate that the 
analysis of Ref. is applicable to the case of radion 
oscillations even though the equation of state is different 
in the case of the radion. However, the coupling of the ra- 
dion to the trace of the matter energy momentum tensor 
opens up a new possibility for achieving tracking at early 
times. The radion could start very close to its effective 
potential S(f> eq i minimum at early times with subdomi- 
nant oscillations (S(p m ax <C 6(f> eq i) induced by the time 
dependent shift of the minimum of the effective poten- 
tial 13.24f) . The radion energy during this early evolution 
can track the matter energy density. At late times when 
Sfimax > 5cj> eq i the oscillations can dominate and lead to 
accelerating expansion for n < 1. This mechanism has 



the possible advantage of avoiding the introduction of a 
tuned potential. Instead only the existence of a stabiliz- 
ing potential minimum is required and a power law be- 
havior around it. The viability of this type of mechanism 
is an interesting problem worth of further investigation. 

A possible direct observational signature of radion os- 
cillations is the existence of a coherent high frequency 
gravitational wave background. Even though the fre- 
quency of these gravitational waves is expected to be of 
order 10 12 Hz and is much beyond the range of current 
experiments, a detailed study of their features could be 
interesting. A study of gravitational wave backgrounds 
from extra dimensions induced by other mechanisms can 
be found in Ref. 0. 

The stability of the field configurations against spatial 
fluctuations is also an interesting issue. Homogeneous 
oscillating scalar fields have been shown to fragment into 
baryons (Affleck-Dine mechanism |35j), Q-Balls |3(j or I- 
Balls |37j depending on their potential and number of 
components. Since the oscillating radion corresponds to 
a real scalar field it may be possible to fragment into 
a generalization of I-Balls (BD I-Balls) with interesting 
cosmological consequences. 
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